Trapped Ion Chain as a Neural Network: Error Resistant Quantum Computation 



O 

o 

(N 
C 



> 

o 
t> 
o 

O 

9 L,: 
-i— > 

c 

=3 

CT 1 - 



13 



Marisa Pons 1 , Veronica Ahufinger 2 , Christof Wunderlich 3 , Anna Sanpera 4 , 
Sibylle Braungardt*, Aditi Sen(De)*, Ujjwal Sen*, and Maciej Lewenstein 5 
1 Departamento de Fisica Aplicada I, Universidad del Pais Vasco, 20600 Eibar, Spain. 
2 ICREA and Grup d 'Optica, Departament de Fisica, 
Universitat Autdnoma de Barcelona, 08193 Bellaterra (Barcelona), Spain. 
3 Fachbereich Physik, Universitat Siegen, 57068 Siegen, Germany. 
4 ICREA and Grup de Fisica Tedrica, Departament de Fisica, 
Universitat Autdnoma de Barcelona, 08193 Bellaterra (Barcelona), Spain. 
5 ICREA and * ICFO-Institut de Ciencies Fotbniques, 08860 Castelldefels (Barcelona), Spain. 

We demonstrate the possibility of realizing a neural network in a chain of trapped ions with 
induced long range interactions. Such models permit one to store information distributed over 
the whole system. The storage capacity of such network, which depends on the phonon spectrum 
of the system, can be controlled by changing the external trapping potential. We analyze the 
implementation of error resistant universal quantum information processing in such systems. 



The Cirac-Zoller proposal of trapped ion computer [l| 
has become one of the paradigmatic models to imple- 
ment a quantum computer Q . Recently, spectacular ex- 
perimental progress in realization of simple algorithms 
and implementation of quantum logic has been achieved 
using a few ionic qubits (c.f. 3]). Although the achieve- 
ment of an all-purpose quantum computer in the near 
future seems difficult, one can be quite optimistic about 
the applications of chains of trapped ions as quantum 
simulators. Recently, it has been shown that long range 
(LR) pairwise interactions between individual spins are 
induced in an ion trap, when applying an additional 
state-dependent force acting on the ions 0, [B| . Also a 
state-dependent optical force can evoke LR couplings and 
has been proposed to simulate spin 1/2 chain systems 0]. 

Here we show that ion spin systems can serve to re- 
alize a neural network (NN) model. NN are aproto- 
type model of parallel distributed memory 0, [J, and 
have been intensively studied by physicists since the fa- 
mous paper by Hopfield [9]. These disordered systems 
with LR interactions, typically present a large number 
of metastable (free) energy minima, like in spin glasses 
Q. These states can be used to store information dis- 
tributed over the whole system. The patterns stored 
have large basins of attraction in the thermodynamical 
sense, so that even fuzzy ones are recognized as perfect 
ones. For this reason, attractive NN's can be used as as- 
sociative memory. At the same time, NN are robust, so 
that destroying even a large part of the network does not 
necessarily diminish its performance. The above listed 
properties make NN's interesting for distributed quantum 
information (QI), where quantum bits do not correspond 
to the internal states (spins) of individual ions, but to 
patterns of the internal states of the whole chain (all-up, 
all-down, half-up-half-down, etc.). These patterns echo 
the lowest energy vibrational modes of the system. The 
sign of the displacement of each ion with respect to its 
equilibrium position in a given mode fixes the up/down 
state of the spin in the associated spin pattern. Some 



approaches to exploit the potential of NN models for QI 
processing have been discussed [l(|, also with respect to 
entanglement generation Here we propose, for the 

first time, a feasible implementation of NN, and the re- 
alization of distributed QI using a chain of trapped ions. 

We first remind the readers the main features of the 
Hopfield model and discuss its similarities with the 
effective Hamiltonian derived in Ref . 0, 0] , that suggest 
the possibility of using a chain of trapped ions as a NN. 
We find the ion-chain storage capacity and its robust- 
ness the most appealing features of NN for distributed 
QI. Thus, the question of ergodicity and, therefore, the 
ability of the system to act as an associative memory 
is not considered here. We show that the storage ca- 
pacity, which is determined by the phonon spectrum of 
the system, can be controlled by modifying the shape of 
the external axial trapping potential and/or by applying 
longitudinal magnetic fields. Although this is a classi- 
cal property of the network, spin ion systems also permit 
to study quantum NN by applying a transverse magnetic 
field or an optical field that effectively simulates it. Here, 
we exploit the storage capacity of the system to perform 
distributed QI i.e. single and two-qubit gates by applying 
appropriate external axial and transverse magnetic fields. 
Transverse magnetic fields should also permit tunneling 
processes between stored patterns and to realize, for ex- 
ample, quantum stimulated annealing [1, 12|. This study 
is beyond the scope of this Letter and will be treated else- 
where. 

Following the models of Hopfield § and Little [l3[ , a 
neuron can be viewed as an Ising spin with two possible 
states: "up" (S = +1) and "down" (S = —1) depending 
on whether the neuron has or has not fired an electro- 
magnetic signal, in a given interval of time [7[. The state 
of the network of N neurons at a certain time is defined 
by the instantaneous configuration of all the spins {Si} 
at this time. The dynamic evolution of these states is de- 
termined by the symmetric interactions among neurons, 
Jij = Jji. Also, full connectivity is assumed, i.e., every 
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neuron can receive an input from any other one, and send 
an output to it. The Hamiltonian reads: 



N N 



(1) 



where h corresponds to an external magnetic field. The 
interactions are determined by the patterns or configura- 
tions of spins to be stored in the network. These patterns 
will be learned if the system is able to accommodate them 
as attractors, implying that a large set of initial configu- 
rations of the network will be driven dynamically to those 
patterns. A possible choice of the interactions is 
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with i j. The p sets of {£f } = ±1 are the patterns to 
be stored. The network will have the capacity of storage 
and retrieval of information, if the dynamical stable con- 
figurations (local minima) reached by the system {Si} 
are correlated with the learned ones {£f }. Although the 
interactions have been constructed to guarantee that cer- 
tain specified patterns are fixed points of the dynamics, 
the non-linearity of the dynamical process may induce 
additional attractors, the so-called spurious states. 

Recently it has been shown that the Hamiltonian de- 
scribing a linear chain of harmonically trapped ions ex- 
posed to a magnetic field gradient [B[ or interacting with 
convenient laser fields [6| can be transformed into an ef- 
fective spin-spin hamiltonian with LR interactions (Jj"), 
mediated by the collective motion of the ions: 
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stants of the chain [14 1 



S^n , where nfj 



with a = x,y,z, label the ions, a are the Pauli 

matrices, F a the force in the a direction experimented 
by the ions, m the ion mass and io a n the angular fre- 
quency of the vibrational mode n. Mf n are the unitary 
matrices that diagonalize the vibrational Hamiltonian: 

are the elastic con- 
The coefficient Mf n gives the 
scaled amplitude of the local oscillations of ion i around 
its equilibrium position, when the collective vibrational 
mode n is excited. Thus, the eigenvectors of M describe 
each ion's contribution to a given vibrational mode, while 
the eigenvalues provide the frequencies, co a ^ n of the col- 
lective modes. 

The external trapping frequencies are chosen such that 
the laser cooled ions crystallize in a linear chain (i.e. 
u x ,i = w Vt i 3> u Zi i) and the external forces act on the 



z— axis (i.e. F x — F v — 0), so that the index a is 
dropped. Henceforth, we consider zero magnetic fields 
Bi = [15| . If we substitute the Pauli matrices in Eq. ((3|) 
by Ising spins S = ±1 (where the effective spin corre- 
sponds to the internal state of the ion), we recover Eq.® 
and the possibility to implement a classical NN with this 
system arises. Nevertheless, there arc some differences 
between both models. First, in the Hopficld model, the 
interactions (Eq.©) are determined by the patterns to 
be stored {£f } = ±1, while in the trapped ion chain, 
the interactions are fixed by the collective modes of the 
system, i.e., the coefficients Mi n that do not necessarily 
equal ±1. Second, in Eq.®, p corresponds to the num- 
ber of patterns to be stored, which in the limit of large N 
(number of spins), is bounded from above by p = 0.14iV 
[7|. In Eq.®, the sum extends over the total number of 
vibrational modes which is larger than the total number 
of stored patterns (spin configurations that the system is 
able to recover). And finally, in the Hopficld model all 
the patterns have the same weight while in the ion chain 
each vibrational mode is weig hted by 1/loI (Eq.®). To 
reproduce as closely as possible a NN behavior, the most 
relevant requirement is the degeneracy of the vibrational 
modes. Moreover, the corresponding patterns must have 
large basins of attraction, i.e. they should correspond 
to sufficiently different spin configurations, so that each 
one is dynamically recovered, even if several spins are 
randomly flipped. 

To check the feasibility of implementing a NN model 
in these ion spin systems, we first find the phonon spec- 
trum using a standard diagonalization procedure, and 
impose the learning rule i.e., we calculate the spin inter- 
actions Jij , mediated by the collective modes of the ions 
(Eq.®). Then, we map a given vibrational mode into 
a spin configuration (initial spin configuration), evaluate 
its energy (Eq.®) and check its dynamical stability un- 
der spin flips using a standard Metropolis algorithm in 
a classical Monte Carlo code. This stability is essential 
for adiabatic QI processing. Explicitly, from the initial 
configuration we randomly flip r spins, and let the sys- 
tem evolve towards equilibrium assuming a noiseless sce- 
nario. If the system recovers the initial configuration, it 
is stable under the flip of r spins. We define the initial 
overlap as rtii — (N — r) /TV. After dynamical evolution, 
the final overlap is given by rrij — (N — s)/N where s 
is the number of spins that differ from the initial con- 
figuration. We repeat this process over M initial con- 
figurations each with r random spin flips. We evaluate 
statistically the final overlap with the initial configura- 
tion as: to/ = (X^^Li m j n j)/M i being nj the number 
of times that the system reaches the configuration with 
777^ . The value of m, for which significant decrease of 
777/ occurs, is a good measure of the size of the pattern's 
basin of attraction. For the harmonic trapping potential, 
the two lowest vibrational modes are the center of mass 
(CM) (all spins parallel, with loi) and the breathing (B) 
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FIG. 1: (a) Ratio between the frequencies of the second and 
first vibrational modes as a function of the exponent of the 
trapping potential for 20 Ca + ions, (b) Final overlap aver- 
aged over 500 initial configurations as a function of the initial 
overlap, for the patterns associated with the two lowest vi- 
brational modes of 40 ions in a potential V = p|a;| ' 5 with 
p = 6.6 x 10~ 20 J/m 1 / 2 . The black squares (triangles) corre- 
spond to the first (second) pattern. 



mode (half up, half down, with ll>2 = V3u>i) [lfj]. We 
find that the pattern associated with the CM mode is 
stable up to the flip of half of the spins, while the one 
associated with the B mode is already unstable under a 
single spin flip. Thus, only the spin configuration associ- 
ated with the CM mode can be recovered (i.e., stored). 
To increase the storage capacity of the network, we con- 
sider here V(x) = p \ x achievable using additional 
control by dc electrodes ([J|(b), see also [17|). We cal- 
culate the ratio o^/^i, as a function of 7, for Ca + ions. 
For N > 20, this ratio depends neither on the number 
of ions, nor on the value of p. For 7 > 1, the ratio is 
rs \/3, and as in the harmonic case, severe limitations 
on the storage capacity appear (see FigQJa)). However, 
for 0.25 < 7 < 0.8, the two lowest modes become nearly 
degenerate. The storage capacity for a system of 40 Ca + 
ions trapped in a potential with 7 = 0.5 is displayed in 
Fig. D^b), where the final overlap is depicted as a func- 
tion of the initial overlap for the patterns associated with 
the two lowest vibrational modes, is close to 1 up to 8 
spin flips, meaning that the system is able to recover four 
patterns (the two associated with the two lowest modes 
plus the two corresponding to a global spin flip). The 
system sometimes reaches a slightly deformed configura- 
tion, which differs only in 1 spin flip from the original 
one (spurious states), making rrif slightly smaller than 1. 
Specifically, the probability of recovering the two modes 
is above 98%, up to 3 initial random spin flips, and above 
97% up to 8 (m, = 0.8). 

Having shown that our system can be robust for clas- 
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FIG. 2: (Color online) Fidelities of the H gate (left) and the 
Bell gate (right), with respect to time (r2 = 0.95ri). The 
classical fidelity bounds are 2/3 and 2/5, respectively (hori- 
zontal dashed lines). The inset in gray shows the fields AX 
and BiX = W 5 BX and B 2 X = 10~ 6 BA with respect to time. 
For adiabaticity, the chosen fields require T > 7x 10 6 ft/A. 
The fidelities are largely independent of the actual dynamics 
of the fields. 



sical information storage, we explore now its capability 
for distributed quantum computing in an error resistant 
way (cf. fill), i.e. robust with respect to the partial 
damage of the system. To this aim, we consider a system 
of 8 spin- 1/2 particles (as in [lj|) in a trapping potential 
with 7 ~ 0.5. The vibrational spectrum is, except for the 
lowest two modes, highly non-degenerate, with rapidly 
increasing eigenvalues. Thus, we consider only the con- 
tributions of the spin configurations associated with the 
three lowest motional modes which, up to a gauge trans- 
formation, correspond to: all-up, half-up-half-down, 2- 
up-4-down-2-up. We encode the information in the first 
two, and consider the third one as noise. The phonon 
mode amplitudes are approximated as M i>m = £- n = ±1 
(exact for periodic boundary conditions). Additionally, 
time dependent "magnetic" fields in the z (B±, B 2 ) and x 
(A) directions are applied leading to the following expres- 
sion for a Quantum Neural Network (QNN) Hamiltonian: 



HQNN(t) = 

+A(S% + S| 

+ s. 



n(SZ + S z 2 + SI + Slf 



z\2 



s-Siy + rz{Si-S z 2 -Sl+Si) 



B 2 {S^ + 5*4 



(5) 



where S° = ofi-i + (T 2ii « = 1, • • • ,4, and r x w r 2 > r 3 > 
0. With A in energy units, all the other parameters in 
Hqnn are dimensionless. B\, B 2 , and A are chosen ini- 
tially (t = to) such that the ground \G(t )}, and the first 
three excited states \Ej(to)) correspond to: | TTtTTTTT)^- 

1 uuuu),, 1 mmu) 2 , and i mmm*. 

To demonstrate universality we focus first on a single 
qubit operation. Identifying |0) = |G(t )), |1) = |-Ea(*o))> 
we consider a single distributed-qubit operation, with a 
generic qubit ao|G(*o)) + di\Ei(to)), evolving adiabati- 
cally under the changes of the magnetic fields, to a final 
state a \G(t = T)) + a x \Ei{t = T)). We choose the time 
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dependence of the fields so that the final state is approx- 
imately the H-rotated state of the input one, where TL 
acts on the logical states as |0) — * |+), |1) — > — |— ), 
where |±) = (|0) ± |l))/\/2. This onc-qubit gate TL is 
reminiscent of the Hadamard gate that takes |0) — > |+), 
|1) — > |— ). The TL operation is achieved by changing 
the positive initial values of B\ and B 2 adiabatically to 
zero, and simultaneously increasing the zero initial A to 
a positive value much larger than the initial B\ , B2 (in- 
set Fig. [2j. Since we deal with superpositions of en- 
ergy eigenstates, we consider in the adiabatic transport 
the dynamical, as well as the Berry phases [2(|. The fi- 
delity of the TL gate is shown in Fig. [^left) as a function 
of time for different noise ratio r 3 /ri. Note that artifi- 
cially increasing the ratio r 3 /ri imitates inaccuracies in 
the trapping potential, disturbance in the motion of the 
ions, as well as noise in the spin (as the phonons are the 
carriers of interaction between the spins). The fidelity is 
quite insensitive to high noise levels, and is larger than 
the classical (measure and prepare) bound of 2/3. 

Let us move now to the two-qubit gates, and treat 
the 4 left spins as one qubit, and 4 right ones as an- 
other, so that: |00) = | TTTTTtTt)*, |oi) = | TTTTUU)*, 
etc. (We have checked that the TL gate fidelity is ro- 
bust in this encoding also.) We demonstrate here, a 
way to implement an entangling universal gate (2ll | act- 
ing as |00) -> (|00) + |11))/V2, |il) - (— 100> + 111))/ y/2, 
|01) -> (|01) + |10))/V2, |10) (-|01) + |10»/V2. We 
denote this gate as Bell gate. We now encode an arbi- 
trary two-qubit state Ylij=o 1 a v l*i)> a oa\G(to)) + 
anl-Bi(to)) + a 01 \E 2 (to)) + a' w \E 3 (t )). The same vari- 
ation of the magnetic fields as in the TL gate leads now 
to the Bell-gate rotated state (Fig. [2] (right)), with fi- 
delity that is noise insensitive and surpasses the classical 
bound of 2/5. Note that in addition to being resistant 
against noise induced by increasing r^/ri, the fidelities 
are also robust against spin-flip errors, as we have en- 
coded the qubit (s) in the two (four) lowest energy levels, 
for the TL (Bell) gate, which we have already shown to be 
metastable against spin- flips. The time scales for which 
both TL and Bell gate fidelities reach maximum values 
are long enough to ensure robust implementation and 
also robustness against errors in time of observation. 

Summarizing, we have shown that spin-ion systems 
can be used to implement NN models. We have calcu- 
lated their storage capacity and robustness against spin 
flips as well as their dependence on the trapping poten- 
tial. Identifying the qubits with configurations of spins 
that echo the lowest vibrational modes of the system, we 
have shown that the system can perform error resistant 
universal distributed QI processing. We have demon- 
strated that by applying adiabatically-varying time de- 
pendent "magnetic" fields, the system realizes single and 
two distributed-qubit operations in a robust way [22l |. 
The scalability issue is like other proposals and experi- 
ments in ion-trap quantum computing [23J , and may po- 



tentially be overcome by connecting mesoscopic clusters 
of trapped ions for instance, by "flying" qubits. 
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